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Under  the  subject  grant  titled  “Analysis  and  Modeling  of  ASREX  Results”  research  was 
conducted  on  analysis  and  modeling  of  low  frequency  sea  surface  scattering  data  obtained  during 
the  Acoustic  Surface  Reverberation  Experiment  (ASREX).  An  important  part  of  the  surface 
scattering  problem  is  the  development  of  accurate  scattering  models  for  the  mean  scattering  level 
as  a  function  of  environmental  conditions.  It  is  equally  important,  however,  to  understand  the 
statistical  distribution  of  the  backscattered  intensity.  The  main  part  of  the  research  addressed  this 
issue  through  modeling  studies  and  through  analysis  of  ASREX  data.  The  final  report  for  this 
research  is  given  below. 

BACKGROUND 

The  Acoustic  Surface  Reverberation  Experiment  was  carried  out  during  the  period  from 
December  1993  to  March  1994  in  the  North  Atlantic.  During  ASREX  Harry  DeFerrari  and  co- 
workers  at  the  Uftiversity  of  Miami  measured  ambient  noise  and  low  frequency  sea  surface 
backscattering  at  short  range.  In  support  of  these  acoustic  measurements,  extensive 
environmental  measurements  were  made  by  other  investigators.  The  development  of  a  better 
understanding  of  low  frequency  sea  surface  scattering  through  the  analysis  of  ASREX  data  was 
the  final  focus  of  the  surface  component  of  the  Acoustic  Reverberation  Special  Research 
Program  (ARSRP). 

For  several  decades  the  Chapman-Harris  empirical  model  [Chapman  and  Harris,  1963]  for 
average  backscattering  strength  as  a  function  of  wind  speed,  frequency,  and  grazing  angle 
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summarized  knowledge  of  low  frequency  sea  surface  backscattering  strength.  Recently,  the 
Critical  Sea  Test  (CST)  program  has  produced  an  improved  empirical  model  for  surface 
backscattering  strength,  the  Ogden-Erskine  model  [Ogden  and  Erskine,  1992],  and  a  large  body 
of  useful  information  on  low  frequency  (100  Hz  -  1  kHz)  surface  scattering.  Important  basic 
research  questions  on  surface  scattering,  however,  still  remain  to  be  answered.  While  it  is 
generally  agreed  that  rough  interface  scattering  is  the  dominant  scattering  mechanism  at  the  low 
end  of  the  frequency  range  (near  100  Hz),  and  that  scattering  from  bubble  clouds  becomes  the 
dominant  mechanism  at  higher  frequencies  and  wind  speeds,  the  environmental  factors  that 
determine  the  spatial  distribution  of  air  void  fraction  near  the  surface,  which  determines  the 
strength  and  statistical  distribution  of  the  scattering,  are  not  well  understood.  The  Ogden-Erskine 
model  is  based  on  wind  speed,  which  is  the  dominant  environmental  factor,  but  clearly  other 
(presently  unknown)  environmental  conditions  are  important  as  well. 

In  addition  to  improving  our  ability  to  predict  scattering  strength,  it  is  also  important  to  improve 
our  understanding  of  (and  develop  models  for)  the  statistical  distribution  of  the  surface  scattered 
returns,  particularly  in  the  upper  end  of  the  distribution,  since  these  returns  can  lead  to  false 
targets.  As  yet,  models  for  the  surface  scatter  statistical  distribution  have  not  been  developed 
from  CST  measurements,  though  studies  by  Gauss  and  by  Huster  show  highlights  with 
lognormal  intensity  distributions  in  some  cases.  In  the  absence  of  other  information,  the  surface 
scattered  field  statistics  would  likely  be  assumed  Gaussian,  which  means  the  scattered  intensity 
distribution  would  be  taken  as  exponential.  However,  lognormal  intensity  distributions,  in 
comparison  to  the  exponential  distribution,  would  lead  to  very  different  probabilities  of 
occurrence  for  intensities  well  above  the  mean. 

ASREX  differed  from  CST  surface  scattering  experiments  in  several  respects.  Probably  the  most 
important  was  that  the  equipment  was  moored  below  the  surface  and  operated  continuously  for 
nearly  three  months  with  source  transmissions  every  12  minutes.  This  mode  of  operation  is  very 
good  for  identifying  environmental  dependencies  and  for  statistical  distribution  studies. 

There  are  important  basic  questions  that  arise  when  the  issue  of  the  backscattered  intensity 
distribution  is  addressed.  When  the  statistics  are  non-Gaussian,  the  issue  cannot  be  simply 
handled  by  measuring  the  intensity  distribution  as  a  function  of  environmental  conditions,  and 
then  using  that  information  directly  in  a  model.  The  fundamental  problem  is  that  the  statistical 
distribution  depends  on  the  scattering  area  (or  equivalently,  on  the  number  of  scatterers)  that 
simultaneously  contributes  to  the  received  intensity.  This  means  that  the  transmitted  pulse  length 
and  the  scattering  geometry  will  enter  into  a  prediction  of  scattered  intensity  statistics.  A 
lognormal  intensity  distribution  may  result  for  a  certain  scattering  patch  size  due  to  intermittency 
in  the  spatial  distribution  of  bubble  clouds,  but  if  this  area  could  be  made  sufficiently  large  to 
include  a  very  large  number  of  such  clouds,  Gaussian  statistics  will  prevail  as  a  result  of  the 
Central  Limit  Theorem.  These  considerations  are  especially  important  for  ASREX,  because  the 
scattering  range  is  relatively  short,  leading  to  a  relatively  small  scattering  patch  size.  In  order  to 
apply  probability  distribution  information  obtained  from  ASREX  to  longer  range  geometries,  a 
methodology  must  be  developed  to  account  for  differences  in  scattering  patch  size. 
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Work  on  modeling  intensity  statistics  was  done  in  two  phases.  In  the  first  phase,  an  analytic 
model  was  developed,  based  on  assuming  that  in  the  scattering  area  there  are  N  (with  N  an 
integer)  independent,  identically  distributed  scatterers.  For  preliminary  comparisons  with 
ASREX  data,  the  individual  scatterers  were  assumed  to  have  a  lognormal  intensity  distribution. 
The  equations  defining  this  model  are  contained  in  the  viewgraph  presentation  given  at  NRL- 
SSC  in  October  1995;  a  copy  is  included  at  the  end  of  this  report.  The  numerical  results  showed, 
as  expected,  that  the  resultant  intensity  distribution  for  N  scatterers  evolves  from  near  lognormal 
to  near  exponential  as  N  varies  from  small  to  large.  In  addition,  the  preliminary  comparisons 
with  ASREX  intensity  distribution  data  showed  that  by  appropriately  choosing  the  model 
parameters  (the  standard  deviation  of  the  individual  scatterer  lognormal  intensity  distribution, 
and  the  number  of  scatterers  N)  good  agreement  between  data  and  model  results  was  obtained. 
(See  “Modeling  Surface  Scatter  Intensity  Statistics”  by  E.I.  Thorsos,  D.B.  Percival,  and  K.M. 
Bader,  presented  at  Surface  SRP  Meeting,  October  17-18, 1995;  included  as  Appendix  A.) 

The  good  agreement  in  these  data-model  comparisons  turned  out  to  be  somewhat  fortuitous, 
since  the  initial  ASREX  intensity  distributions  were  produced  by  grouping  intensities  over  the 
entire  experiment  in  particular  grazing  angle  bins.  Thus,  the  environmental  conditions  varied 
considerably  over  the  time  period  used  to  construct  the  distributions.  In  addition,  the  data 
processing  algorithms  introduced  some  intensity  averaging  for  each  transmission,  which  would 
also  affect  the  intensity  statistics. 

For  phase  two  of  this  work,  improvements  were  made  in  the  data  analysis  and  the  model.  Neil 
Williams  at  the  University  of  Miami  reprocessed  800  Hz  surface  scattering  data  specifically  for 
intensity  distribution  analysis,  and  all  averaging  was  removed  from  the  processing  algorithms. 
The  raw  data  was  then  supplied  to  APL-UW  for  further  analysis.  We  carefully  removed  all 
obviously  spurious  data  points  and  grouped  the  data  in  both  wind  speed  and  grazing  angle  bins. 
Unfortunately,  questions  still  remained  at  the  end  of  this  effort  on  the  quality  of  the  ASREX  data 
set.  It  was  apparent  that  significant  overall  level  uncertainties  existed  in  the  data  set.  While  it  is 
possible  that  the  intensity  distributions  relative  to  the  means  are  still  valid,  uncertainty  in  this 
regard  led  to  the  decision  not  to  publish  the  results  of  this  analysis. 

The  model  was  improved  in  the  phase  two  analysis  by  removing  the  rather  artificial  assumption 
that  the  number  of  scatterers  in  the  scattering  area  was  an  integer  N.  More  realistically,  the 
average  number  of  scatterers  <N>,  not  necessarily  an  integer,  is  specified  together  with  a  model 
for  the  pdf  of  N.  If  one  makes  the  simplest  assumption  that  the  scatterers  are  spatially 
uncorrelated,  then  it  follows  that  the  number  of  scatterers  in  the  scattering  area  should  be  Poisson 
distributed.  Using  this  assumption  an  analytical  model  for  the  intensity  distribution  has  been 
developed,  and  a  Monte  Carlo  implementation  of  the  model  was  also  been  made.  The  two 
models  give  excellent  agreement  (see  viewgraphs  discussed  below). 

Results  of  the  phase  two  analysis  were  presented  at  the  Acoustical  Society  Meeting  in  December 
1996,  and  copies  of  the  viewgraphs  are  included  at  the  end  of  this  report.  [Thorsos,  Percival,  and 
Bader,  1996;  see  Appendix  B.]  For  this  intensity  statistics  analysis,  only  data  points  with  SNR  > 
15  dB  were  used.  No  intensity  distribution  variation  was  noted  over  the  grazing  angle  range  from 


10  deg  to  30  deg,  and  therefore  the  intensity  distributions  over  this  grazing  angle  range  were 
combined.  The  scattering  area  on  the  surface  ranged  from  20,000  to  40,000  square  meters. 

The  ASA  viewgraphs  summarize  the  results.  At  low  wind  speed  (up  to  5  m/s)  the  intensity 
distributions  were  found  to  be  fit  very  well  with  an  exponential  distribution.  This  means  the  field 
statistics  are  Gaussian.  At  these  low  wind  speeds,  the  scattering  should  be  from  the  rough  sea 
surface  and  scattering  from  bubbles  should  not  be  important.  The  exponential  distribution  for  sea 
surface  scattering  implies  that  the  rough  sea  surface  has  a  large  number  of  scatterers  within  the 
scattering  area.  At  high  wind  speeds  the  intensity  distribution  deviates  from  exponential  with  a 
higher  upper  tail.  This  is  the  regime  of  bubble  cloud  scattering. 

The  model  described  above  was  used  to  fit  the  high  sea  state  distributions.  This  can  be  done 
assuming  the  scattering  area  contains  a  set  of  identical  scatterers  that  scatter  with  a  lognormal 
intensity  distribution,  and  the  number  of  scattering  is  Poisson  distributed.  Good  agreement  is 
found  with  eight  scatterers  on  average  for  the  ASREX  geometry,  and  for  the  scatterer  intensity 
distribution  to  have  a  standard  deviation  of  5.0  -5.5  dB.  To  fully  test  such  a  model,  data  sets  are 
needed  for  a  range  of  scattering  areas. 

REFERENCES 

Chapman,  R.  P.  and  J.  H.  Harris,  Surface  backscattering  strengths  measured  with  explosive 
sound  sources,  J.  Acoust.  Soc.  Am.  34, 1592-1597, 1962. 

Ogden,  P.  M.  and  F.  T.  Erskine,  An  empirical  prediction  algorithm  for  low-frequency  acoustic 
surface  scattering  strengths,  Naval  Research  Laboratory  technical  report  NRL/FR/5 160-92-9377, 
April  28  1992. 

Thorsos,  E.  I.,  D,  B.  Percival,  and  K.  M.  Bader,  “Modeling  low-frequency  surface  scatter 
intensity  statistics,”  J.  Acoust.  Soc.  Am.  100, 2799  (1996) 


Eric  I.  Thorsos 


cc:  ONR  Administrative  Grants  Officer  (June  Hawley),  1  copy  memo,  copy  of 

Form  298,  copy  of  DTIC  Card 

Director,  Naval  Research  Laboratory,  Code  2627, 1  copy  memo,  enclosures,  and 
Form  298 

Defense  Technical  Information  Center,  2  copies  of  memo,  enclosures,  298  Form,  and 
DTIC  50  Card 

Grants  &  Contracts,  University  of  Washington  (Sinh  Simmons),  1  copy  memo 


Modeling  Surface  Scatter  Intensity  Statistics 

Appendix  A 


Eric  l  Thorsos 
Don  B .  Percival 
Kate  M.  Bader 


Applied  Physics  Laboratory 
College  of  Ocean  and  Fishery  Sciences 
University  of  Washington 
Seattle,  WA 


Surface  SRP  Meeting 
October  17-18, 1995 
NRL-SSC 


Approved  for  public  rlease;  distribution  is  unlimited. 


Introduction 

•  Scattering  strength  average  scattered  intensity 

•  Need  model  for  intensity  pdf 

•  Cannot  simply  measure  intensity  pdf  for  direct  use 
in  model 

-  Intensity  pdf  depends  on  the  scattering  area 
(i.e.,  on  the  number  of  scatterers) 


Modeling  issue:  Predict  intensity  pdf  for  large  scattering 
area  when  pdf  is  measured  for  small  scattering  area 
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Gaussian  Statistics 


•  Notation: 

Complex  field :  X  +  iY 
Magnitude :  R  =  a/x2  +  Y2 
Intensity :  I  =  R2  =  X2  +  Y2 

•  When  large  number  of  scatterers  contribute  to  field 

-  X  and  Y  are  Gaussian  distributed 

-  R  is  Rayleigh  distributed 

-  I  is  exponentially  distributed 

•  In  this  case,  pdfs  are  independent  of  geometry 


For  ASREX,  I  is  not  exponential 

—  Closer  to  lognormal 
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Intensity  Level  PDF 


Exponential  pdf: 

P,(I)  =  i  e-''» 


•  Lognormal  pdf: 

where  r|  is  mean  and  o  is  standard  deviation  of 
normal  distribution  for  /nl 


•  Intensity  level: 

L  =  10  log!  o  (1/(1)) 

Intensity  level  pdf  pdf  of  L 

[for  convenience,  (I)  =  1] 
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Model  for  Intensity  PDF 

•  Need  pdf  for  total  intensity 
I  =  X2+Y2 
where 

X  +  iY  =  I(Xn  +  iYn) 

n=1 


•  Assume 

-  Xn  +  iYn  is  independent  of  Xm  +  iYm  (n  ^  m) 

-  Contribution  from  each  scatterer  has  lognormal 
intensity  distribution 


•  Need  joint  distribution  PX  Y(X,  Y)  for  single  scatterer 

Let  X  =  R  cos  $ 

Y  =  R  sin  <]) 

•  Assume 

PfuM  =  pr(r)p^(^)  =  ^:pr(r) 
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Model  for  Intensity  PDF  (continued) 

•  Joint  characteristic  function:  1  scatterer 

OO  OO 

♦x,y(ui-U2)=  (dX  (dYPx,Y(X,Y)ei(ulX+U2Y) 

— oo  — OO 

•  Joint  characteristic  function:  N  scatterers 

<I)X,Y,n(u1>u2)  =  [‘t'XjY  (U1  >  u2  )] 

•  Joint  pdf  for  total  field 

OO  OO 

Px,y,n(X,Y)  =  —5-2  ( dU!  } du2  <t>x,Y,N(ui,u2) e“i(uiX+U2Y) 

(2n) 

\  /  — oo  — oo 

total  intensity  pdf 
average  intensity  pdf 
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Model  for  Intensity  PDF  (continued) 

•  Joint  characteristic  function:  1  scatterer 

OO 

^X,y(u)=  JdlPi(I)  J0  (VI  u) 
o 

•  Joint  characteristic  function:  N  scatterers 

^X.Y.nC11)  =  [^X.yC11)] 

•  Average  intensity  pdf 

Kl  00 

Pi,n(I)  =  5  |duu <I>X,Y,n(u)  Jo (VNl  u) 

^  0 
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Remarks  on  Intensity  Statistics  Modeling 

•  Model:  N  independent  scatterers  in  scattering  area 

•  Need  intensity  pdf  for  single  scatterer  as  function  of 
environmental  conditions 

•  Need  number  of  effective  scatterers  per  unit  area  as 
function  of  environmental  conditions 

•  Should  consider  N  as  random  variable 

-  Need  pdf  for  N  per  unit  area 
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Modeling  Low  Frequency 
Surface  Scatter  Intensity  Statistics 
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Introduction 

•  Scattering  strength  average  scattered  intensity 

•  Require  other  statistics,  e.g.,  intensity  pdf 

•  If  acoustic  field  statistics  are  Gaussian,  intensity  will 
have  exponential  pdf 

•  CST  results  (Gauss,  Huster)  at  low  frequency 
(100  Hz  - 1000  Hz)  indicate  intensity  pdf  has  higher 
tail  than  exponential 

•  Cannot  simply  measure  intensity  pdf  for  direct  use  in 
model 

-  Intensity  pdf  depends  on  the  scattering  area,  i.e., 
on  the  number  of  scatterers 

Topics  to  be  discussed: 

•  Measurements  of  intensity  pdf  —  ASREX 

•  Modeling  the  intensity  pdf 

-  Dependence  on  scattering  area 
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Acoustic  Surface  Reverberation  Experiment 

(ASREX) 

University  of  Miami 
Harry  A .  Deferrari  (PA. ) 

Neil  J.  Williams 

Date:  December  93  -  March  94  (84  days) 
Location:  34°  N,  70°  W 
Frequencies:  100,200,400,800  Hz 
Transmissions  every  12  minutes 


Initial  data  analyzed  for  intensity  statistics 

-  Frequency:  800  Hz 

-  Time  period:  first  27  days 

-  Grazing  angles:  10°-30° 

-  Pulse  type:  13-digit  Barker  code 

-  Effective  pulse  length:  20  ms 
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ASREX  Acoustic  Mooring 
(University  of  Miami) 
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Histogram  of  Intensity  Fluctuations 

•  Removal  of  deterministic  trends 

-  Intensities  grouped  into  one-degree  bins 

-  Local  intensity  average  obtained  using  2-hour 
running  means 

-  Intensity  fluctuations  relative  to  local  mean  yield 
histogram 

•  Histograms  conditioned  on 

-  Wind  speed 

-  SNR  >  15  dB 

•  Scattering  area  -  2-4  x  104  m2 
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Probability  Density  Probability  Density 


ASREX  Intensity  PDF 
800  Hz,  10-30  Degrees,  Days  353-380 
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Modeling  Surface  Scatter  Intensity  PDF 

Assumptions: 

•  Effective  discrete  scatterers 

-  Scattered  field  a  random  variable 

•  Scattered  intensities  identically  distributed  for  all 
scatterers 

-  Bubble  cloud  scattering 

-  Use  lognormal  distribution 

-  Need  to  specify  standard  deviation 

•  Scattered  phase  uniformly  distributed 

•  Scatterers  randomly  distributed  spatially 

-  Implies  number  of  scatterers  in  a  given  area 
is  Poisson  distributed 

-  Need  to  specify  average  number  of  scatterers 
per  unit  area 
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ASREX  Intensity  PDF 
800  Hz,  10-30  Degrees,  Days  353-380 
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Summary 

•  Conditions:  f  =  800  Hz 

Grazing  angies  10°-30° 

Scattering  area  -  2-4  x  104  m2 

•  At  low  sea  states  (U  <  5  m/s)  intensity  distribution  is 
exponential 

-  Region  of  interface  scattering 

•  At  higher  sea  states  intensity  distribution  deviates  from 

exponential  with  higher  upper  tail 

-  Region  of  bubble  cloud  scattering 

•  Model  can  replicate  ASREX  intensity  distributions,  but 

-  Model  is  underdetermined 

-  Need  to  compare  with  data  for  smaller  scattering 
areas 
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